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Generation of spherically symmetric metrics in / (R) gravity 

Z. Amirabi0 M. Halilsoylll and S. Habib Mazharimousavlil 
Department of Physics, Eastern Mediterranean University, Gazimagusa, Turkey. 

In D—dimensional spherically symmetric / (i?) gravity there are three unknown functions to be 
determined from the fourth order differential equations. It is shown that the system remarkably 
integrates to relate two functions through the third one to provide reduction to second order equa¬ 
tions accompanied with a large class of potential solutions. The third function which acts as the 
generator of the process is F {R) = ■ We recall that our generating function has been employed 

as a scalar field with an accompanying self-interacting potential previously which is entirely different 
from our approach. Reduction of / {R) theory into system of equations seems to be efficient enough 
to generate a solution corresponding to each generating function. As particular examples, besides 
known ones, we obtain new black hole solutions in any dimension D. We further extend our analysis 
to cover non-zero energy-momentum tensors. Global monopole and Maxwell sources are given as 
examples. 


I. INTRODUCTION 

/ (R) gravity is one of the modified theories of Ein¬ 
stein’s general relativity that attracted much attention 
in recent times [H-I^ . In @ / (i?) = R + aR^ with a > 0 
has been introduced as the model of inflated universe 
while / (R) = R — a/R”' (a > 0, n > 0) was consid¬ 
ered as a candidate for the dark energy model (iSl. 
This model, however, is not a viable model for dark en¬ 
ergy and instead / (R) = R — aR^ with a > 0 and 
0 < n < 1 emerged as alternative which has been pro¬ 
posed i n I12I.I13II . Later on more viable models were stud¬ 
ied in [l4l-[^ . A detailed review of these models are 
given in [l9j| fOther review papers are given in [20l - [^ l. 
Some recent works on solutions in / (i?) gravity are given 
in [ 2344 ^ . For / (i?) = i?, in D = 4 dimensional space- 
time it coincides with the standard general relativity but 
otherwise it stands with an action that is arbitrarily de¬ 
pendent on the Ricci scalar. Finding exact solutions in 
this theory with fourth order derivatives of the metric 
tensor is both important and challenging. Apart from 
exact analytic solutions there are / [R) models that can 
only be expressed implicitly in non-polynomial expres¬ 
sions. Each particular model has advantages / disadvan¬ 
tages as far as experimental tests are concerned [43 - l53 |. 
There are even models that lack the Einstein’s A—gravity 
limit. Among other expectations the UV / IR behaviors 
at near / far distances, quantum renormalizability with 
power counting of the counter terms are prominent. At 
any cost preferring to abide by the classical regime we 
confine ourselves first to sourceless (vacuum) / (i?) mod¬ 
els that admits exact integrals. In the last section we 
extend our discussion to cover external sources such as 
global monopole [sj and electromagnetic field. Let us 
add that the equivalence of / (R) gravity to Brans-Dicke 
(BD) theory {uj = 0) with a potential has also been high- 
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lighted extensively in the past as a transition between 
Jordan and Einstein frames. In this approach the ex¬ 
act solutions can be generated by adopting scalar field 
ansatzes which in general brings into the Lagrangian in¬ 
tricate potentials. Our method will be confined entirely 
to the Jordan frame without reference to the BD field or 
any scalar potential. 

The vacuum of / (i?) gravity is known to carry its own 
curvature sources. By vacuum in this theory it is meant 
the absence of an external energy-momentum tensor T^j/, 
of any physical source . Carames and de Mello in lat¬ 
ter work have considered the spherically symmetric vac¬ 
uum solutions of / (i?) gravity in higher dimensions. We 
shall rederive most of their results anew together with 
some additional extensions which we are interested in to 
develop further. Addition of external ^ 0, no doubt 
makes the problem technically more complicated, but fol¬ 
lowing the lesson learned from the vacuum / empty so¬ 
lutions of / {R) = R theory of gravity we will attempt 
to derive the most general equations and in some specific 
cases the solutions as well. The D dimensional spheri¬ 
cally symmetric line element that we shall consider will 
be 

ds^ =—A{r) di^-\ - \-—dr‘^-\-r'^dUl \,_2 (1) 

B{r) 

in which A (r) and B (r) are metric functions to be deter¬ 
mined while dn ^_2 represents the (D — 2) —dimensional 
unit spherical line element. In particular integrals we 
have the restrictive case A(r) = B (r) included, but more 
general cases with A{r) ^ B (r) must be interesting as 
well. Beside the metric functions A (r) and B (r) we 
shall employ a third function denoted by F [R) = 
which characterizes the type of the / [R) gravity. Let 
us add that not in all cases of / (i?) models the explicit 
form of / [R) can be expressed analytically in terms of 
the variable R, the Ricci scalar. Instead, it involves a 
transcendental part that can’t be inverted in the form 
of r = r (i?), these are hybrid forms. We add that 
even these hybrid forms don’t prevent us from calculating 
^ > 0 and 4^ > 0 which are crucial terms to deter- 
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mine the absence of ghosts and thermodynamic stability, 
respectively. 

In brief what has been achieved in this paper is to 
show that the metric functions A (r) and B (r) are related 
through an integral expression for the function F {R) (or 
F (r)). This amounts to the fact that once F (r) is given 
it acts as a generator to generate a new set of {A (r ), 
B (r)) pair. The function A{r) is expressed in terms of 
B (r) and F (r) and the remaining equation is reduced 
into a master equation satisfied by B (r). Once we give 
an ansatz for F (r) our master equation can be integrated 
in principle to obtain B (r). In this manner we can ob¬ 
tain an infinite class of metrics in / (i?) gravity generated 
from an infinite set of F (r). No doubt the dimensionality 
of spacetime F>(= d+l) also plays role in the derivation. 
In particular, we present examples of new black hole so¬ 
lutions in D > 3, by the method described above. We 
wish to add also that in the reduction process the system 
of differential equations in / [R) gravity reduce naturally 
from the fourth order to the second order. 

The paper is organized as follows. In Section II we red¬ 
erive the / (i?) held equations in ZI—dimensions which is 
comparable in some sense with [s^ . A number of exam¬ 
ples to justify the effectiveness of our method are given. 
Generalization to ^ 0 is analyzed in Section III. We 
end our discussion with Conclusion in Section IV. 


II. THE FIELD EQUATIONS IN 
D-DIMENSIONS 

The dimensional vacuum / (i?) —gravity is repre¬ 
sented by the action 

in which / (i?) is a function of Ricci scalar R and D > 3. 
Variation of the action / with respect to provides the 
held equations (in metric formalism) 

FB!'^ - IfS; - + S-^DF = 0 (3) 

in which F = -^ and □ is the covariant Laplacian. The 
general spherically symmetric line element is given by (1) 
and the held equations (3) are explicitly given by 


in which 


^ =-TT 
* 4A 


B'A' + 2BA” - 


BA'"^ 2D2BA' 
+ —^- 


Rl = -^[ b'A' + 2BA" - 
4A 


A 

BA'^ 2D2B'A 

- A -^- 

A r 


( 7 ) 

( 8 ) 


and 


R2=- 


1 


2r^A 


{rBA' F rAB' + 2 D 3 A {B - 1)), (9) 


for 1 < * < I? 2 . Herein and in the rest of the paper, 
Dk = D—k, a prime stands for the derivative with respect 
to r and F = ^. A particular combination of these three 
equations leads to two equations which are independent 
of /. The hrst equation (4) may be written as 


2ABrF" -h H {AB' - BA') = 0 
which integrates to 

F' 


where 


A = BF[‘^ exp ( —2Di / 'j^dr 


H = rF' F D 2 F ^ 0. 


( 10 ) 

( 11 ) 

( 12 ) 


The second equation (5) upon considering (11) becomes 
independent of A too. The closed form of the second 
equation reduces to a linear equation for B (r) , which is 
given by 


2D-i 

B" F PB' F 2QB F^=0, 


where 


H 


H 2 
rF r 


and 


Q = 


rF'" {H^ - FirFF') F" F-H 


H 


(13) 

(14) 

(15) 


FiZ2 r^F 

Finally, the explicit form of / in terms of r is given by 

2B {2rF' F FD 3 ) F" 


f = 


H 


2rF' [rB' F D 3 B) - 2 (rB' F D 3 {B - 1)) F 


(16) 


To complete our analysis we give the explicit form of the 
Ricci scalar in terms of r which reads 




, / BA'F' D 2 BF' 


and 


( 0 , 


^ ^ _ D 2 { 2 tB' FD 3 {B- 1)) 

r 2 

2r^BF'" ^ r'^F" {3r'^B'F' F FD 2 { 3 rB' F2 BDi)) 
Fr2 F2r2 ’ 

(17) 

In summary, the only equation to be solved is Eq. (13) 
which is second order and linear for B. Therefore the 
procedure is reduced to set a F (r) - which eventually 
represents the form of / (F) - and solve the only equation 
(13) to find B (r) and consequently A (r). 
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A. Applications of the method 


2 . -F (r) = 1 + ar, with a =constant 


Before we give certain applications for our formalism 
we would like to compare our approach with the work 
of Carames and Bezerra de Mello [55|. The main differ¬ 
ence can be seen from the fact that in there are two 
generating functions (so to say) which are F and Y. In 
other words Eqs. (16) and (17) of [HI are coupled and 
one must consider them together to find a solution to 
the field equations. In our formalism we have only one 
generating function which is F and Eq. (13) is the only 
equation to be solved. In the following three cases we 
shall show that for simple cases {F = 1 and F = 1 -|- ar) 
our results overlap with [ssj but for more complicated 
case {F = ar°-) our solution is the general one while the 
solution given in [HI is a restricted one (look at Eq. (46) 
in [ 5 ^ and (35) and (36) in this work). 


1. F{r) = l 


We start with the simplest case, with F = 1 or equiv¬ 
alently FI = D 2 . Definition of F = ^ implies 

f = R-2A (18) 

in which — 2A is an integration constant to be interpreted 
as the cosmological constant. The main equations (11) 
and (13) admit 


B (r) = A (r) 


l + ^ + C2r^, n>3 
Ci+C 2 r^, D = 3 


(19) 


for integration constants Ci and C 2 and consequently 

/(r) = -2C2Fi, (20) 

R = -C 2 DD 1 


with 


As of Sec. 3.2 of Ref. [S^ we consider F" (r) = 0 
in our general formalism but instead of going through a 
general F—dimensional solution, we investigate the cases 
in closed form for / (R). From Eq. (10) with F" = 0 one 
obtains A = B (up to a constant which one can set it 
unity via a redefinition of time). For an arbitrary D the 
solution for B may not be possible in a closed form but 
for specific dimensions we may find. 

a. D = 3 In three dimensional spacetime the solu¬ 
tion for B becomes 


B = A = C 2 r^ + Cl 


and 



Cir^a^ In 



(23) 


/ 


— 4(72 + 4Cia^ In 



2 Cia(l + 2 ar) 

r(l + ar) ^ ’ 


where Ci and C 2 are integration constants with the cur¬ 
vature scalar 


R = - 6 C 2 + 6 Cia^ In 



Cia (2 + 9ar + 6 a^r^) 
r (1 -|- ar)^ 


(25) 

The hybrid relation between / (r (R ), R) and R is given 
by 


3 3r (1 -I- ar) 


(26) 


Note also that from the foregoing expressions we can 
identify the cosmological constant by (72 = — A. Further¬ 
more, setting a = 0 , one recovers the previous example 
with Cl scaled which suggests that it is related to the 
mass of the central object. 

b. D = 4: In four dimensions the solution is given 
by 


A — ——D 2 DiC 2 - 


( 21 ) 


Finally the solution becomes 


1 _ 2M _ 2A 2 n t 
r TTTTDj n„ n, ^ i U > o 


B{r)=A{r) = { ^ (22) 


in which we set 


_2M ]J <1 

(7i = <! Ds ^ 

^ ^ -M D=3 


-I- — (sar — 2 — -|- 60 ;^ In M -|- — ) ) (27) 


with 


/ = -6(72 - 36(7ia'' In 1 -b — 

ar 

6 aCi [-1 + 6 a'^r'^ + 3ar) 2ar + 1 

-\-1-T— (28) 

r^ [1 + ar) r^ 


where M is the ADM mass. The solution for F > 3 
is Schwarzschild de / anti de-Sitter black hole solution 
and for F = 3 it is the BTZ black hole. Let’s add that 
F (r) = ^ =constant, does not change the nature of the 
solution for the metric with constant scalar curvature. 
We note that our results in this section expectedly is the 
same as section 3.1 in [H| . 


and 


F =-12(72 - 72(7ia3 In 1 


1 

ar 


6a(7i (—1 + Qa'^r'^ + 6 ar) (1 -b 2Qfr) 
(1 -b ar)^ 


( 29 ) 
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We see clearly the role oi C 2 = — ^ and we wish to pro¬ 
ceed with Cl = 0 which implies 

B = ^ = + l + (30) 

so that 

f = R + 2aVR - 4A - 2A. (31) 

This is a black hole solution with a singularity at r = 0 
such that 

R=L+4.A. (32) 

If we set A = 0, with a < 0 we may introduce a horizon 
for the solution located at 


We also note that although a and a are two arbitrary 
constants a must satisfy a ^ —D 2 , 1 ± y/Di- This is re¬ 
markable to observe that in Sec. 3.3.2 of the same 
ansatz for F has been considered but the solutions to 
the field equations (see Eq. (46)-(49) of [H) are not the 
same as what we found here is more general. As a matter 
of fact our solutions (35)-(38) with C2 = 0 reduce to their 
solutions. This shows that reducing the field equations 
into a master equation with a single generating function 
makes advances in finding exact solutions in / (i?) grav¬ 
ity. 

In H = 3 dimensions the solution becomes rather spe¬ 
cific since the last term vanishes for all values of a. The 
functions then read as 

2a{a-l) 

A = r^^B, (39) 

-2 q 2 2(l + 2a-a.^'j 

B = Cir^^+C2r ^ , (40) 


while for a > 0 the solution possesses a naked singularity 
at r = 0. To complete our investigation let’s determine 
the absence of ghosts and thermodynamic stability of the 
explicit f(R) found in (31). We see that ^ = \ + —^== 

and = — 2(-j; both conditions can not 

be satisfied simultaneously. 


/ 


4aC2 f2a^ — a — l) ao-a) 

_ 1 _ L'p a + l 

n -f 1 


„ 2 (2a-b 1) (a - 3) C2 2a(i-a) 

R = — -- - -r “+i . 


(41) 


(42) 


This is nothing but the solution found by Zhang, Liu and 
Li in [5^ with their parameters 


Our next example is a power-law form for F (r ), i.e., 
F = (34) 

with constants a and a which upon (11) yields 

2o(a-l) 

A = r “+-“2 B. (35) 

Substituting into (13) one finds 


2a2-6a+6 + (2a-5)D + D2 2(i32+2a-a2) 

B = CiT “+^2 -b C2r “+^2 

113(04-112)^ 

(2a2 - 6a -b 6 -b (2a - 5) H -b (II2 -b 2a - o^) 

(36) 


in which Ci and C 2 are the integration constants. Fol¬ 
lowing A and B we also find 


2aC'2lli (a — 1) (II2 + 2a) r “+°2 2 aaDiD^r°‘ ^ 

a -b LI2 D 2 -b 2a — a^ 

(37) 

and 


R = 


C 2 D 1 {D - a) {D 2 + 2a) 


(a -b ZI2) 


r 


2a.(,o.-l) 
0. + D-2 


aPiDs (g - 2) 

(ZI2 + 2a - a2) ^ > 


P = 


2 a (1 — a) 
0- -|- 1 


(43) 


and 


kL = — 


2 (2a -b 1) (a - 3) Ci 
a -b 1 


(44) 


Note that 


/ ~ R^^ (45) 

which yields f ^ R? for the specific choice a = 

For H > 4 one may set C2 = 0 and therefore 


_ 2a^-6a+6+(2a-5)P + P^ 

B = Cir “+°2 

113(04-712)^ 

(2a2 - 6a -b 6 -b (2a - 5) H -b {D 2 + 2a- a?) 

(46) 


with an analytic relation for / (i?) given by 

/ = aR}-^ (47) 

in which d is a constant which can be set to unity (by a 
fine choice of a). In the case of a = 0 the theory gives 
i?—gravity. The solution is a black hole with singularity 
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at r = 0. Let us also add that with Ci = 0 and a = 1 the 
solution reduces to 


ds^ 


-^dt^ + \dr‘^ 


^dVl 


D 2 


in which 


e = 


Ri 

D 2 


(48) 


(49) 


with 


4- A new black hole solution in D = 3 

In [13 where / (i?) = in three dimensional space- 
time with d = const, has been studied, the solution does 
not cover the case d = —^ which makes / (i?) = VR- This 
can be seen from the Eq. (12) of and Eq. (45) (note 
that 2 (i~-a) ~ 5 answer) and in that paper both 

do not cover the case f {R) = VR- However, in what 
follows we wish to show that the solution / (R) = \fR 
can easily be obtained in three dimensions. To do so let’s 
consider 


and 



f = '/R. 


(50) 


F (r) = jdr exp (ar) 


with a and (3 two real constants. This yields 


(51) 


A=— + C2r\ 

r 


(58) 


(59) 


This represents a global monopolc-type solution with a 
deficit angle. Another interesting setting is for a = —2 
and C 2 = 0 which upon a proper choice of a one finds 


B = 



exp (—2ar), 


(60) 


f{R)=R^ (52) 

and the solution becomes 


/ = 2a/3 ( — -I- 4 C2r^ ) exp {—ar) 


(61) 


DDq+26 

B (r) = Cir 


D3DI 

{DDq -|- 26) Diq 


(53) 


while 


A{r) =r^B{r). 


(54) 


Clearly £> = 10 and D = A are excluded. D = 4 is not 
allowed directly from (10) where H = 0 with F = -^ 
demanding AB = 0 which is not acceptable. Eor D = 10 
the particular solution can be obtained as 


B = C 2 


Qj.6 



(55) 


and 


and 


R = Aa 


Rl 

4r2 


C2r 



exp (—2ar). 


(62) 


Note that, Ci and C 2 are two integration constants such 
that C 2 effectively plays the role of a cosmological con¬ 
stant. The solution is a black hole with a horizon located 

at Th = the condition that > 0. Also 

from R we see that the solution is singular if and only if 
Cl jtz 0. In the sequel we are interested in (72 = 0 which 
makes the solution to be rather simple but singular. Ac¬ 
cordingly the form of / and R are given by 

f o exp(-ar) 

f = 2a(3Ci -, (63) 

r 


/ 


4a (9(72 - 7 - 21 In r) 


with the Ricci scalar 

„ 98-b 168 In r - 72(72 


(56) 


(57) 


The metric function A{r) follows accordingly from (11) 
which is given in (54). 

We comment that / = aR^~^ does not satisfy ^ > 0 

and > 0 simultaneously unless we set a to be negative 
(Note that d = 1 is needed to have the Einstein gravity 
recovered.). For instance with d = 1 and a = —2 which 
implies f = R^, and both conditions are satisfied. 


R = aCi 


exp {—ar) 


(64) 


which upon tuning the free parameter fd by 4a,5^(7i = 1 
the form of / becomes 


f = RR- 


(65) 


Here a and (7i are positive constants and the line element 
finally reads 


ds^ 


^- dr^ + r^dO^ 

r To exp (—2ar) 


( 66 ) 


in which tq = is also a positive constant. 
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5. A general class of solutions in 2 + 1 — dimensions 

In three dimensional spacetime in addition to what 
we found by now, we wish to show that there exists an 
important class of solutions yet to be discovered. This 
specific class, however is a characteristic feature of only 
three dimensions. To see this solution let’s set ZI = 3 in 
(13) which yields 


The solution is a black hole solution with a < 0 and 
therefore one may write 

(l - ^) 

in which the horizon is shown as . A change of variables 
of the form t = V^T, r = -j= reduces (76) to 


r^F {rA' - 2A) F" 

- r [{rA" - A') F + {rA' - 2A) F'] {rF' + F) = 0 (67) 

where we substituted B (r) by A{r) using (11) i.e.. 


Ae^ t-f'+f 
{rF' + Ff 


( 68 ) 


Eq. (67) possesses a trivial solution for A (r) irrespective 
of the form of F (r) which is given by 

A {r) = Cgr'^ (69) 


in which Co is an integration constant. In this situation 
the field equations are all satisfied provided B (r) and 
A{r) satisfy the condition (68) i.e.. 


ds^ = - (l - dT^ + + \p^dn^ (77) 

which is the Schwarzschild black hole with a deficit angle 
caused by a cosmic string. 


III. GENERALIZATION TO / {R) GRAVITY 
COUPLED TO MATTER SOURCES 

In this Section we extend our vacuum analysis to the 
presence of matter coupled with gravity. Therefore the 
action becomes 

1 = 1 d^xV^ (^ 8 ) 


C'or^e'^'^ rF^+F'^’’ 
{rF' + Ff 


(70) 


One can see easily that with F = ^ =constant we obtain 


in which Cm is the matter Lagrangian density. The field 
equations become 

FR"^-\f5';-^"^^FF5';UF = 8iTGT"^ (79) 


R C’o 2 

B = ■ 

(71) 

The solution given by (69), (71) and 


f{R) = p + ^R 

(72) 


with rj = const.^ constitutes a particular class. Note that 
since the choice of F (r) in (70) is arbitrary this can be 
used to generate an infinite class of solutions. This will 
not be searched any further here. 


6. 3 + 1 —dimensional black hole solution in 

f {R) = R + 2a\/H gravity 


Previously we found an exact solution for the model of 
gravity in the form 

f{R) = R + 2aVR (73) 

with a 7^ 0. The solution to the field equations is given 
by 

A(r) = B(r) = i + ^ (74) 

with 

i?=l. (75) 


with = diag {—p,p, q, q) the energy momentum tensor 
of the matter source. The line element is going to be a 
spherically symmetric as (1) and without going through 
the details of the field equations, we give the changes in 
the field equations (11) and (13). The corresponding field 
equation to Eq. (11) reads 

/l = Ba^exp(-2 /(k^-k^)*) (80) 


while the main equation (13) i.e. the master equation for 
B (r) takes the form 


B'‘ 


P 


r{p + p) 
HB 


B' + 2QB + 2S- 


2 r^ {p + p) 


in which P, Q and H are given in (14), (15) and 
respectively and 


■ = 0 

(81) 

( 12 ), 


Do , {p-q) {DiF' -rF"){p + p) 

= ^ + — --■ (82) 


Note that, unlike the source-free case, here the master B 
equation is not a linear equation. It can also be observed 
that the particular choice of p -I- p = 0 removes the non¬ 
linearity in the B equation (81). Further choice oi p = q 
leaves us with the same equation for B as in the sourceless 
case. Yet with q ^ 0, the source shows itself in the d' 
function as given in the sequel. 
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The closed form of / is given by 


/ 


- 2 ( 7 - 


r'^A 




where 


and 

(83) f:^R + 2C2 + ^^+0{a^) (91) 

with 


= 2r^ABF" + [BrA' + {rB' + 2 BD 3 ) A] rF' 

- F [BrA'+ {rB'+ 2 {B - 1) D 3 ) A] . (84) 

The foregoing expressions are not impressive much unless 
we provide concrete examples. This is our aim in the 
section that follows. 


A. Applications 


1. 3 + 1 — dimensional global monopole eoupled to 

f — R + 2ay/R modified gravity 


One immediate application of Eq. (81) is the extension 
of the / = i? + 2a\fR in 3 + 1—dimensional vacuum to 
the gravity coupled to the global monopole whose 
energy momentum at very large distance is given by (we 
assume A = B vn the line element) 

r; = d*a(7.^[l,l,0,0] (85) 


in which 77 represents the global monopole charge. The 
solution for the metric is given by 


A = B = 


+ 1 ) 

2 >ar 


( 86 ) 


with 



(87) 


and 


F (r) = 1 + ar. (88) 

Again we must have a ^ 0 and in the limit 77 = 0 one 
recovers the vacuum solution. 


S. 2 + 1 — dimensional Maxwell electric field coupled to f{R) 

gravity 

Let’s consider now Maxwell electric field to be coupled 
to / {R) gravity in 2 + 1—dimensions. Using the line ele¬ 
ment (1) and the standard Maxwell Lagrangian together 
with F = 1 + ar one finds 

= diag.^ (89) 

in which Q is the electric charge. The solution with small 
a up to first order simply reads 

A = B~2Q‘^ (2ar-l)lnr 

+ 2a{Q^ + M)r-M + -h O (a^) ( 90 ) 


R ~ -GCa + ^ 

- ^ ( 4 Q 2 + M + 2Q‘^\nr) +0 (a^) . (92) 

Clearly by setting a = 0 one recovers the charged BTZ 
black hole with cosmological constant A = C 2 . We add 
that the solution when the energy momentum tensor of 
the matter source is of the form of a fluid, one should 
consider it as an interior solution. Therefore one has to 
make sure that the Israel junction conditions are satisfied 
at the interface between exterior and interior solutions 
[s^. In two examples we studied here, the energy mo¬ 
mentum tensors are long range fields which allow us to 
consider our solutions to be exterior. 


IV. CONCLUSION 

The integrability of / (i?) vacuum gravity with spher¬ 
ical / circular symmetry is reduced first to a set of mas¬ 
ter relations, i.e., Eq. (11) and a master equation for 
B (r ), i.e., Eq. (13). Given any ansatz generating func¬ 
tion F = -^ in terms of the coordinate r, our method 
generates a solution pair of A(r), B (r) and a solution 
for the function / (i?). Most of our solutions encountered 
are of hybrid nature, that is / (i?) can’t be expressed ex¬ 
plicitly in terms of R. The power-law form for instance, 
of the form f (R) ~ i?^, with a rational number fc, is 
obtained easily in our method. Some of the solutions 
presented as applications are already known. Yet, new 
and rare type of solutions can also be obtained easily. 
In the second part of the paper, we extend the integra¬ 
bility of vacuum case to the non-vacuum / (R) theories. 
For this case we found also a master relation, i.e., Eq. 

(80) and a master equation, i.e., Eq. (81). Consider¬ 
ing F (R) to be the generating function for our formal¬ 
ism and Tff = diag [—p,p, q, g] to be our energy momen¬ 
tum tensor one finds a solution to the master equation 

(81) . We presented two examples. In the first example 
we set F = 1 -I- ar with a global monopole coupled to the 
gravity in 3 -I- 1—dimensions. The second example con¬ 
siders the Maxwell electric field coupled to the gravity 
in 2 -|- 1—dimensions with the same generating function. 
In this case we found the solutions approximately for 
small a. Let’s also add that among the few examples we 
have studied we found some closed form of / (F) such as 
f{R)=R + 2ay/R - 2A - 2A and / (F) = dF^-t . We 
found that the first one can not satisfy the conditions for 
the absence of ghosts and thermodynamic stability simul¬ 
taneously while the second one for specific a does satisfy. 
Finally it will be in order to state that our method of re¬ 
duction for the spherically symmetric / (F) gravity has 





a large scope as far as solutions are concerned. Physi¬ 
cal implications of the solutions obtained such as dark 
matter connection is not considered in the present study. 


Extension of our formalism to problems with different 
symmetries, such as stationary axial symmetry requires 
a separate investigation. 
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